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models. We consider overall scale factor perturbation of Einstein spaces as well as more 
general anisotropic deformations of three-spheres. At leading order in a', renormalization- 
group equations turn out to be Ricci flows. In the three-sphere background, the latter 
is the Halphen system, which is exactly solvable in terms of modular forms. We also 
analyze time-dependent deformations of these systems supplemented with an extra time 
coordinate and time-dependent dilaton. In some regimes time evolution is identified with 
renormalization-group flow and time coordinate can appear as Liouville field. The resulting 
space-time interpretation is that of a homogeneous isotropic Friedmann-Robertson- Walker 
universe in axion-dilaton cosmology. We find as general behaviour the superposition of a 
big-bang (polynomial) expansion with a finite number of oscillations at early times. Any 
initial anisotropy disappears during the evolution. 
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1. Introduction 

Finding string backgrounds with cosmological interpretation is a challenging problem that 
has led to many developments, either in genuine string theory or in string-inspired mod- 
els. One possible way to investigate string cosmology is by promoting a Euclidean three- 
dimensional static string solution to a Minkowskian four-dimensional time-dependent one. 
To some extent, this is borrowed from general relativity, where e.g. a three-sphere can 
be promoted to a Friedmann-Robertson-Walker (frw) universe. The latter is a four- 
dimensional space-time which has spherical spatial section of time-dependent radius. How- 
ever, string theory is more constraint than general relativity and such a promotion is a del- 
icate issue, which does not admit a systematic treatment: the string equations of motion 
involve many fields (metric, dilaton and antisymmetric tensors) and the four-dimensional 
universe must necessarily be accompanied by an internal, compact space that saturates the 
central charge without spoiling the perturbative (in a') nature of the solution. It is not 
surprising therefore that the emergence of e.g. de Sitter solution in the context of strings is 
still unclear, whereas this universe is the archetype of a cosmological background in general 
relativity. 
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One can be more concrete and give a specific example to illustrate the above difficulties. 
A three-sphere background appears in string theory as the target space of a Wess-Zumino- 
Witten (wzw) model. This sigma-model is an exact rational conformal field theory. Com- 
pared to the principal model, it contains a Wess-Zumino term, which is topological and 
gives rise to a Dirac monopole quantization (put differently, the Kalb-Ramond three-form 
is classified according to the third homotopy group of the group manifold, which in the 
compact case is always tts^G) = Z). The radius is therefore quantized and this creates an 
obstruction when trying to let it vary continuously while keeping conformal invariance. 

A possible way out would be to abandon conformal invariance. The perturbed sigma- 
model at hand is no longer an infra-red fixed point and various parameters (such as the 
radius of S* 3 in the above example) now are running with the renormalization-group (rg) 
energy scale. It is tempting therefore, to interpret the RG scale as a cosmological time 
and the various (new) fixed points as steady or final states of the cosmological evolution. 
Such identification is however questionable because time is expected to appear as a genuine 
dimension of the target space of a critical string and nothing a priori guarantees that its 
evolution is the RG-flow of a non-critical string with one dimension less. Indeed, string 
backgrounds satisfy second-order equations while RG-flow equations are first-order. Hence, 
there is no reason to trust the above approach generically, although the dissipative nature 
of the dilaton may help in some circumstances to transform second-order into first-order 
equations. In some other examples, off criticality can generate a Liouville mode that 
may eventually play the role of the time direction, putting thereby on a firmer ground 
the cosmological interpretation of the RG evolution (see for instance Refs. [1, 2, 3, 4, 
5]). Despite these reservations, the method remains useful in order to acquire a wider 
perspective of the RG landscape around a given spatial string solution 1 . 

Instead of abandoning criticality, one can try to keep it or restore it. Starting from a 
conformal sigma-model with spatial target space, we add an extra time direction (which 
requires adjusting the internal manifold) and let some parameters explicitly depend on it. 
Time evolution of these parameters is dictated by the critical string equations of motion. 
This procedure may go under the name of dynamical promotion of the RG-flow although, 
depending on the parameter and the regime, this evolution may or may not be the one 
obtained in the absence of time, by following the RG-flow. An extreme situation correspond 
to the case where the parameter is a coupling constant associated with a truly marginal 
deformation of the original conformal sigma-model. Such a coupling does not run under 
the RG whereas it can acquire in general a non-trivial time dependence. 

In the present note, we analyze both the RG-flow properties and its dynamical promo- 
tion in the following framework. We start with a conformal sigma-model with space-like 
(i-dimensional target space. We assume the metric be Einstein with constant curvature, 
while criticality is achieved with a Kalb-Ramond field and no dilaton. The off-critical 
excursion is realized by switching on an overall scale factor c for the metric that does not 
affect the antisymmetric tensor. We perform the RG analysis and show, at lowest order in 
a' , that no other fixed point exists away from c = 1. 

*A renormalization-group approach to cosmology has been put forward by other authors [6]. That 
perspective is however slightly different from ours. 
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We further investigate the case of the three-sphere. There, in a mini-superspace ap- 
proximation, we consider a three-parameter perturbation written in terms of left-invariant 
currents describing a squashed three-sphere sustained by a Kalb-Ramond field. The RG- 
flow equations turn out to be equivalent to the Halphen system that we can explicitly solve 
in terms of Eisenstein series. The modular properties of the latter allow to demonstrate 
that, no matter the initial conditions, the system flows towards the round sphere. Again 
no other fixed point exists away from the SU(2) wzw model. 

It is remarkable that in all cases at hand - global scale factor or anisotropic perturba- 
tions, the RG-flow is in fact a Ricci flow driven by a connection with torsion, due to the 
Kalb-Rammond background. Ricci flows have attracted much attention recently [7], not 
only in mathematics but also in the physics literature, in relation to various issues of string 
theory such as tachyon condensation or time dependence (see e.g. [8, 9, 10]). 

We investigate time-dependent solutions by considering a family of backgrounds ex- 
pressed as the warped product of the original (i- dimensional manifold with an extra time 
direction, while keeping the B field unperturbed 2 . The one-loop Weyl-invariance equations 
are satisfied for the full (d + l)-dimensional space-time provided a background dilaton is 
also switched on in the time direction. In the case of the three-sphere, we generalize the 
time warping by implementing time dependence in all three parameters that were studied 
in the RG approach - anisotropic warping. 

The kind of solutions we obtain when only the overall scale factor is made time- 
dependent (simple warping) is by construction of the Friedmann-Robertson- Walker type. 
The system exhibits a friction behaviour. Whenever the latter dominates in the solution, 
we recover the RG-flow of the d- dimensional sigma-model studied previously. In this regime, 
the off-critical string generates a Liouville mode which becomes the "missing time" with a 
background charge that is identified with the friction coefficient. Warping coincides then 
with Liouville dressing. As a consequence, the dilaton of the full (d+ l)-dimensional target 
space is linear in this regime. 

In general, the solution turns out to describe a universe undergoing a power-law expan- 
sion (big-bang solution) superimposed with a damped oscillation governed by the dilaton 
and the axion. This kind of cosmologies was considered e.g. in [11, 12]. 

The case of more anisotropic warping that is studied for the sphere can be summarized 
as follows: the equations of motion are found to be a generalization of the Halphen system 
with a friction term. The solution exhibits the oscillating behaviour met previously. In 
particular, in the linearized limit, around the wzw point, all three parameters oscillate in 
phase and the solution remains qualitatively the same. 

The paper is organized as follows. The RG-flow analysis is performed in Sec. 2 first 
for the general case with perturbation of the overall factor and then for the three-sphere 
with more parameters for RG motion. In Sec. 3 we introduce a time coordinate, in the 
manner described previously and take the various coupling constants perturbed so far to 
be time-dependent. This allows us to write down the equations of motion which are a 

2 This ensures that in the case of compact group manifolds, the quantization condition appearing in 
the WZW point is satisfied despite of the warping and is no longer an obstruction for finding new critical 
solutions. 
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generalization of the equations for the RG-flow and take the form of the motion of a point 
particle in presence of a friction term. The system is non linear and we were not able 
to find a closed solution, but both numerical and asymptotic analyzes are possible. They 
show that the universe we describe grows as a superposition of a big-bang expansion and a 
damped oscillation, and eventually converges to a linear dilaton solution [13, 14]. Section 
4 is actually devoted to the proper cosmological analysis of the solution, performed by 
moving to the Einstein frame. A summary follows in Sec. 5. 

2. The renormalization-group-flow approach 

2.1 The behaviour of the overall scale in general situations 

We would like to study the RG-flow for the coupling of the metric in a system without 
dilaton, i.e. for the sigma model 

£(c) = \ (c 9iiV + B^) dx»8x u , (2.1) 

where g^ v is a d- dimensional Euclidean metric and c and A are constants. The antisymmet- 
ric tensor is tuned so that for c = 1 the model is conformal 3 and for simplicity we require 
g^y to be Einstein with constant curvature (that is always true e.g. for group manifolds, 
but this case is not exhaustive). Geometrically this means that H = dB is the torsion that 
parallelizes the Riemannian metric g^ u , and no dilaton is required for conformal invariance. 

Two-dimensional sigma models are perturbatively renormalizable. We will analyze 
the evolution of a certain class of deformations with special emphasis to their target-space 
interpretation, in the mini-superspace approximation. We will perform this analysis at one 
loop although it is in principle tractable at higher order. 

Following [15, 16, 17, 18] we define the connection 

r ""- = {r p }-H < 2 - 2 ' 

and the corresponding Riemann tensor 

R ~^ispa = ( 1 - ^) R>l u P a- ( 2 - 3 ) 
In a dimensional-regularization scheme the one-loop counterterm reads: 

SC m = £-R- l ,8X"BX- (2.4) 

where R is the Ricci scalar for g^ v . In particular, if we start with an SU{2) wzw model 
at level k, we have d = 3 and R = 6/k. 

J One should not confuse c with the central charge. 
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From the above expressions, we can determine the behaviour of c with respect to the 
scale [x, captured in the corresponding beta-function, /3(c) = dc /diog^. We find: 



which exhibits an infra-red fixed point at c = 1 that we already know. It is possible to go 
beyond one-loop and resum higher-order corrections. In the particular case of wzw models, 
where the target space is a group manifold, this amounts to a finite renormalization of the 
radius, or equivalently, a shift of the level of the affine algebra k — > k + g*, where g* is the 
dual Coxeter number of the algebra. 

It is interesting to remark that given the form of the one-loop counterterm, the evo- 
lution equations take the form of a Ricci flow where the effect of the B field is taken into 
account by the connection r~. In fact the equations above are equivalent to 



dcg^ = J_ R - 
d log jjL 2-ir 



R~nv. (2.7) 



As usual one should pay attention to the direction of the flow and it is hence useful to 
define an RG time variable t = — log n so that g^ v (t) describes the evolution of the system 
going towards the infra-red. 

In order to compare the result with what we will find in the following, we introduce 

ait)= l 2 \ogcit). (2.8) 

Then the evolution of a(t) going toward the infra-red gives: 

with V(a) given in Eq. (3.14). Equation (2.9) admits the implicit solution 



2vrd 



(eMi) _ argtanh (e 2a ®)) + est. (2.10) 



To summarize, in d dimensions and for constant-curvature metrics, the running of the 
one-parameter deformation that maintains the normalization for the topological term is a 
Ricci flow, and it does not produce any new non-trivial fixed points. 

2.2 Generic perturbations of S 3 : the Halphen system 

The case of wzw is of special interest. The underlying algebraic structure allows to switch 
on more general perturbations, while it remains possible to analyze their simultaneous 
evolution under the RG flow. Although the issue of metric perturbation has attracted 
much attention in the past (e.g. for 0(3) or 0(4) sigma models [19, 20]), we will here 
investigate a class of perturbations, for which the full analysis of the RG flow can also be 
carried out in the presence of a Kalb-Ramond field. 
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Let us for the moment consider a WZW model on a generic compact semi-simple group 
G. Its action is as in Eq. (2.1) with c = 1. The corresponding algebra q has generators 



{ T a } and structure constants / . The standard Killing form is 

ds 2 = 5 a(3 J a ® J 73 , 



(2.11) 



where J a are the left currents J a = y/ktr(T a g 1 dg) for g E G. We will consider a deformed 
wzw model with metric 



ds 2 = g aP J a ®J p = Y^^) ja ® JQ 



(2.12) 



(/x is the RG-scale) and the standard B field, which is fixed, being a topological term. It is 
convenient to go to the vielbein of the currents so that in particular the metric reads: 



/ 7i (A*) 



(9a/3) 



72 {lA 



\ 



\ 



7»(m)/ 



(2.13) 



where n = dimG and 7c(a*) arbitrary positive functions. 
The RG-flow equations are here 



1 



dg Q/ 3 _ 

d log [i 2-ir 



1 



ad 



2vr 



R- 



a(3, 



(2.14) 



where R a p are the components of the Ricci tensor, H = dB and = H ai $Hp t& '. We 
observe that the RG-flow equations for the perturbation pattern at hand are again governed 
by a Ricci flow with torsion, as in Sec. 2.1. 

Although the full analysis is tractable for any compact group G, we will here focus on 
the SU(2), where, as we will see below, the flow equations can be solved explicitly. The 
metric (2.13) has now only three entries: 7i(a*),72(m)> 73 (/•*)• Using (A. 9), we obtain the 
following Ricci tensor: 



(Raff) 



/ 7? ~(72 -73) 
7273 

7|-(7l~73) 2 
7173 

V 



73 -(71 -72) 



(2.15) 



71 72 



while the Kalb-Ramond term reads: 



{Hip) = 2 



1 

7273 



1 

7371 



1 

71 72 



(2.16) 



4 We will absorb the level k of the algebra in the currents. Therefore, in terms of the structure constants, 
H — fai3-fJ a A / A J 1 , where the first indices in / are lowered with the Killing metric (H is the standard 
field for the unperturbed wzw model). Its square is contracted with the new metric though. Explicitly, 

Ha/3 = ■ 
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We introduce as before an RG-time pointing towards the infra-red, 

1 



dt 



-dlog/x, 



(2.17) 



2^71(^)72(^)73 (a*) 

where we have also reabsorbed the product of the three 7 a 's. These are positive functions 
and such a rescaling is therefore harmless since it does not spoil the monotonic evolution of 
t with n (t increases from minus infinity to plus infinity when flowing to the infra-red from 
infinite to zero fx). Putting everything together one obtains the following RG equations: 



,71 
71 

72 



(72-73) -7i + l) 



< 2^ = (73 -71) - 7 2 2 + l, 

72 

2^ = (7i-72) 2 -7l + l, 
L 73 



(2.18) 



where dot denotes the derivative with respect to t. 

The evolution of Ricci flows on general homogeneous or locally homogeneous spaces 
has been studied from various perspectives (see e.g. [7, 21]). In the absence of torsion, the 
last constant term in (2.18) is missing and the flow converges towards the round sphere, of 
vanishing radius though. The presence of torsion does not alter this behaviour but affects 
the radius of the sphere which stabilizes to \fk because all 7 a 's now converge to one 5 . This 
non-trivial infra-red fixed point corresponds to the SU(2)k wzw model. 

The above results on the convergence of the flow are based on asymptotic analysis. 
However, as already advertised, the Ricci-flow equations (2.18) can be solved explicitly in 
the case at hand. Indeed, setting 



t = log(T + T ), 



(2.19) 



which amounts in identifying T £ [0, +00 [ with a monotonic function of the inverse energy 
scale 1 /n, and 

7371 



7273 



T + T ' 
equations (2.18) are recast as: 



ft 2 = 



T + T ! 



7172 
T + T '' 



(2.20) 





= n 2 n 3 


-fii (n 2 + n 3 ), 


< tl 2 


= n 3 n! 


-n 2 (n 3 + fii), 


n 3 




- n 3 («! + n 2 ) , 



(2.21) 



where the dot now denotes the derivative with respect to the new time 6 T. The arbitrary 
constant Tq does not play any significant role. 



J Note that the differential systems under consideration are parabolic, hence not invariant under time 
reversal. 

6 In the case of a system without torsion one would obtain the same equations by defining Q.\ — 72 73 
and cyclic and keeping the same time variable t. 
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This is the celebrated Halphen system that was studied in the 19th century, also called 
Darboux-Halphen because the equations were written by Darboux on the analysis of triply 
orthogonal surfaces [22] and solved three years later by Halphen [23, 24]. It has appeared 
since then in several instances in physics 7 , as e.g. in the search for self-dual or anti- 
self-dual reductions in four-dimensional Euclidean gravity (Bianchi IX class with SU(2) 
isometry) [25] or in the scattering of £[/ (2)- Yang-Mills monopoles [26]. The function 
Y = -2 (fii + n 2 + O3) satisfies 



Y = 2YY -3Y 2 , 



(2.22) 



known as the Chazy equation [27, 28]. Note that integrable ordinary differential equations 
like Chazy's or Halphen's turn out to be systematically related to self-dual Yang-Mills 
reductions. 

The Halphen and the Chazy equations possess remarkable properties. If LO a (z) provide 
a solution to the Halphen system for generic z € C (the dot stands then for d /dz), then so 
do 

^ {z) ~ ^) + (2 ' 23) 

with (jj) £ PSL(2, C). Similarly, if y(z) is a solution of the Chazy equation, so is 

1 /„.,^\ 6c 



[cz + a) 

A class of solutions are expressed as [29] 



( az±b\ _ 
\ cz+d ) 



cz + d 



(2.24) 



w a (z) = - ^^log^ a (2:), 



(2.25) 



where E a form a triplet of modular forms of weight two for T(2) C PSL(2,Z) (see e.g. 
[30, 31]), transforming thereby as 



1 

z 





(° 1 ° ^ 






-»2 2 10 




E 3 J 


v° -v 





(2.26) 



and 



z+l : 









(E, 


E 2 


-1 


."11 


E 2 






10 1/ ' 


\E 3 



(2.27) 



The class of solutions under consideration does not capture situations where some w a 's are 
equal. In our physical set-up, this would correspond to deformations of the three-sphere 
that preserve some extra isometries, whereas all-different w a 's restrict the isometry group 
to SU(2), strictly. We will discuss this possibility towards the end of the present section, 
but the following should already be stressed: Eqs. (2.21) imply that if Q = fig then 



7 The connection to Ricci flows was also noticed by K. Sfetsos in unpublished work. 
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£l a = £lp and therefore £l a and £lp never cross unless they are equal at any time. For the 
moment, we will focus on the situation where Q,\ / / ^3- We can express the E a as 

d^/dz d^/dz d-^/dz 

^ 1 = ^T' ^ 2= A^T' ^ 3= A(A^l)' (2 ' 28) 
and recast the Halphen system in terms of the Schwartz equation for A: 

*y*. »^y*"V_ 1 1 WW' (2 . 29) 



dA/dz 2 v dA /d* / 2 va 2 (a - 1) 2 A(A-i);v^y 

The elliptic modular function 

A = I (2.30) 

is a celebrated solution for this equation (see App. B for a reminder on Jacobi theta 
functions). This provides a particular solution of the Halphen system in terms of the three 
r(2) weight-two Eisenstein series. Correspondingly, 

y = 12-^- log r/, (2.31) 

where r\ is the Dedekind function. 

For our present purposes, we must focus on real solutions ft Q for real time T (z = iT). 
These are obtained as 

Q a (T) = iiu a (iT) = -~ log E a (iT) (2.32) 

and 

Y(T) = iy{iT) = A log El E 2 E 3 . (2.33) 

The modular properties of the functions under consideration set stringent constraints be- 
tween the asymptotics of the solution and its initial conditions: large-T and small-T regimes 
are related by T <-> i/t. Assuming 0° = £l a (0) finite and using (2.26) and (2.32), the 
asymptotic behaviour (T — > 00) is found to be universally 

Q a = — + subleading. (2.34) 

This provides an elegant proof of the universality of generic SU(2) Ricci flows in the 
presence of torsion towards the corresponding WZW infra-red fixed point. 

One can further assume fi^ finite and positive. This is natural for describing the 
initial deformation of a three-sphere and sufficient to show that all Q a remain positive at 
any later time. Indeed, suppose that < < < ^3 an d that fti has reached at 
time ti the value Q\ = 0, while > 0. From Eqs. (2.21) we conclude that at time 

ti, £l\ = ^3 = — ^2^3 < and Cl\ = S7g ^3 > 0. This latter inequality implies that fii 
vanishes at t\ while it is increasing, passing therefore from negative to positive values. This 
could only happen if 0^ were negative, which contradicts the original assumption. However, 
if indeed ft? < and ftg, fi§ > 0, there is a time t\ where fii becomes positive and remains 
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Figure 1: A generic solution of the Halphcn system for positive < fiSj < dis- 



positive together with 0,2 and until they reach the asymptotic region, where they all 
satisfy (2.34). The generic behaviour of a positive-initial-value solution is given in Fig. 1. 

Although the Halphen system can be solved, its deep nature makes it difficult to 
establish the correspondence between a given set of initial conditions and the modular 
forms E a (see Eqs. (2.25) or (2.32)) necessary to provide the actual solution. Moreover, 
it is important to stress that solutions exist, for which the initial conditions are not all 
positive and not even finite. Hence, for these solutions, (2.34) does not hold. This happens 
e.g. for the particular solution given in Eqs. (2.28) and (2.30), where T = is a simple pole 
with positive (unit) residue for ^2,^3 and a double pole with negative residue (—V 2 ) f° r 
ill. As a consequence, is negative and increases, at large T, exponentially towards zero, 
while is positive and decreases exponentially towards zero; ^2 is positive and decreases 
exponentially towards n/2. Solutions with negative JTs were considered in Refs. [26, 32] for 
the description of the configuration manifold of two SU(2) monopoles. 

The existence of poles is generic for all solutions of the Halphen system. For solutions 
corresponding to a set of finite and different initial values, these poles are pushed behind 
T = 0. This is related to the following general property: Halphen's solutions possess a 
natural movable boundary [29]. They exist in a domain of C, where they are holomorphic 
and single valued, and this domain has a boundary that contains a dense set of essential 
singularities. The precise location of this boundary depends on the initial conditions. 

We would like finally to mention that the Halphen system admits also solutions which 
are not based on elliptic functions. These solutions have power-like dependence on time: 



They are actually the most general solutions with Q2 = ^3 and describe axisymmetric 
deformations of the three-sphere 8 , namely deformations preserving an SU(2) x U(l) . The 

8 Notice that the equations at hand being first order, setting $1° = ^3 guarantees that Q.2 = at any 
subsequent time. 




(2.35) 
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constants (A, C) are arbitrary and determined by the initial conditions. This class is closed 
under PSL(2,C) transformations. Indeed, using (2.23), we learn that under T — > 1 /t, 
(A,C) -► (i/a,-C/a 2 ), while for T — > T + 1, (A, C) -► (A + 1,(7). Regularity requires 
= 1 /a > 0, which ensures that the pole of VL a is located at T < 0. Furthermore, 
the asymptotic behaviour is again universal as in Eq. (2.34). No new fixed point appears 
therefore in this case either. The RG flow of perturbations of this type, namely with 
7i 7^ 72 = 73) wa s analyzed in [20], without torsion though. When translated into the 
language of self-dual four-dimensional Euclidean metrics with SU{2) isometry, solutions 
(2.35) correspond to the general Taub-NUT family, including Eguchi-Hanson metrics [25]. 

The case of VL\ = VL2 = ^3 corresponds to fully isotropic deformations studied in the 
previous section. In this case, the solution is unique: VL a = 1 /t+a. 

To summarize and conclude the present analysis, solutions of the Halphen system fall 
in three classes: Hi / fl2 7^ ^3> ^1 / ^2 = ^3 an d fii = Q2 = ^3. The corresponding 
three-manifolds, target spaces of the sigma-model, are homogeneous with isometry group 
SU(2), 577(2) x £7(1) and 5*7(2) x SU(2), respectively. This captures all possible isometry 
groups for a general three-sphere, the latter, most symmetric case corresponding to the 
usual round sphere. The Ricci flow describing the renormalization of the sigma-model 
leads unavoidably to the round sphere, which is therefore the unique perturbative infra-red 
fixed point, found at one loop. 



3. Space-time solutions 

3.1 Equations of motion for the dilaton and the scale factor 

As advertised in the introduction, our aim is now to describe a generalization of the previous 
construction obtained by introducing an extra time dimension and treating the coupling 
as a time-dependent field. This clearly bears many resemblances with the usual Liouville 
dressing of [1, 33, 34, 35]. In other words, we would like to write down the Weyl-invariance 
equations for the following sigma model: 

s= l -Jd 2 z [-dtBt + (c(t) 9 ^ + b^) dx»Bx v + m(t)] , (3.1) 

where, g and B are background fields solving the d-dimensional equations of motion. Mod- 
ulo the time-dependent dilaton, the system in Eq. (2.1) appears thus as a description of a 
constant-time slice. 

Let us start by rewriting the d + 1 dimensional metric in the form of a Weyl rescaling, 
9u, = e 2CT « ( ° ) = e^g MN , (3.2) 

where c(t) = e 2a ^\ as before. This means in particular that the Ricci tensor can be written 
as 

-f^MN = Rmn ~ 9mnK l — (d — 1) K MN , (3-3) 
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where K MN is defined as 

K M " = -8 M ag^d L a + 5 NL (d M d L a - T p ML d P a) + ^g LP d L ad P a5 M \ (3.4) 

After some algebra one finds 

Rtt = -d(a(t) + a 2 (t)) , (3.6a) 

Rtp = 0, (3.6b) 

Rpv = R»v + g^u (d& 2 (t) + a(t)) . (3.6c) 

The (d + l)-dimensional metric is not Einstein, so we expect the presence of a time- 
dependent dilaton. The equations of motion read: 

'fo = -ln* + v M $v M $ - ^#MNPtf MNP = ^=^, 

/3 G = 7?mn - + 2V M V N $ = 0, (3.7) 



4"MN 

2 



/?B = — h^ M H MNP + 9 M $iJ M NP = 0. 



Several remarks are in order here. 

• As already pointed out, the action in Eq. (3.1) does not describe by itself a critical 
string theory. An internal conformal field theory of central charge c\ must be super- 
imposed to the above sigma-model. Requiring that we are in the perturbative regime 
with respect to a' , we expect: (i) the model in Eq. (3.1) to have a central charge of 
the order of d + 1 and (ii) the internal CFT with c\ k, 25 — d to have a geometrical 
interpretation with a (25 — d)-dimensional compact target space (we refer here to the 
bosonic string). 

• If the scale factor were chosen to be kept constant and equal to its critical value c = 1, 
the natural solution to the above equations would be the linear dilaton $ = ±qt. This 
solution is actually an exact CFT with central charge 1 — 3q 2 < 1 (the minus sign 
is due to the time-like direction). The central charge Cd of the d-dimensional sigma 
model (decoupled in this case) depends on the geometry (the curvature) of the target 
space. For some spherical-like geometry, one expects < d and therefore the internal 
CFT must have c\ > 25 — d. This imposes a locally hyperbolic internal target space 
that can be e.g. a quotient of a hyperbolic plane H 2 5-d by some discrete subgroup. 

The /3$ = const, equation can be seen as a consequence of the Bianchi identities and 
the (3b = and fie = equations. Explicitly it reads: 



The Pb = condition is satisfied by our choice of Kalb-Ramond field (see Sec. 2.1) and 
a time-dependent dilaton $ = 3>(i). So, let us consider the remaining (3q = equation. 
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First of all we separate the time component (all terms retain a block-matrix structure): 

(3.9) 



R tt + 2d t 8Mt) = -d (a(t) + a 2 (t)) + 2$(t) = 0, 

(1 - e~ 4CT «) + (dd 2 {t) + a(t) - 2&(t)*(tj) = 0, 



where we used the equations of motion for the a = conformal system 

1 



Riiv - (3.10) 



Taking the trace with g MN one obtains the system: 
d (er(t) + a 2 (t)) - 2*(t) = 0, 

ife-M*) (1 - e- 4(J W) + d (da 2 (t) + er(t) - 2&(t)$(tj) = 0. 



(3.11) 



Introducing 

Q(t) = -m+p(t), (3.12) 



the equations become: 



Q(t) = -la 2 (t), 

a(t) = -f e- 2(T W (1 - e~M*)) - 2a(t)Q(t). 



(3.13) 



The second equation has precisely the structure of the equation of motion of a point 
particle in a potential 

V(a) = ^e- 6CT (1 - 3e 4CT ) + est., (3.14) 

in the presence of a time-dependent friction coefficient Q(t). In the Q — > oo limit we recover 
the same equation as in Eq. (2.9) with the same potential V(a) when identifying the energy 
scale jJL for the running, off-shell {i.e. non-conformal) , d- dimensional system with the time 
direction according to 

lirt 

log^ = -^-. (3.15) 

This was to be expected since what we find is the usual link between the RG-flow equations 
and the conformal condition typical of a Liouville dressing (see e.g [1, 2, 33, 35]). This 
identification holds, however, only in the large-Q regime. 

A useful reformulation of the equations is obtained by multiplying the second equation 
by a and introducing the energy E = a 2 / 2: 

f«JM = -aw, (3 . 16) 

|£(t) = -FMt))+!Q 2 (t), 
which in turn allows us to write down an equation for Q as a function of a: 



Ta Q^) = -d\j— d J - - (3-17) 
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(a) Small Q(0) 



(b) Large Q(0) 



Figure 2: Typical behaviour for a(t) in the non-linear system. For (a) small and (b) large (positive) 
initial values of Q(t) (numerical integration with (a) (Q(0), <r(0), <r(0)) = (0.5,0.1,0.1)) and (b) 



3.2 Linearization 

The system (3.13) can be solved numerically and typical results for large Q(0) and small 
Q{0) are shown in Fig. 2. Since we describe a dissipative motion, the system relaxes to 
equilibrium with or without fluctuations. 

The solutions we found are in principle only valid at first order in a' so one may wonder 
about their stability. Actually, it turns out that if the spatial part is supersymmetric (which 
is the case if the model in Eq. (2.1) is of the Wess-Zumino-Witten type), then so is the 
spacetime in Eq. (3.2). This is immediate once one notices that the Weyl tensor for g is the 
same as for the Cartesian product with line element ds 2 = — dt 2 + g^dx^dx" and therefore 
vanishes identically. 

Around the conformal point. A further step can be made via linearization. To study 
the behaviour in a neighbourhood of the a = conformal point we introduce 



(Q(0),ct(0),ct(0)) = (10,0.1,0.1)). 



E(t) 



ait) 



(3.18) 



and derive the first-order system 



Q(t) 

a(t) 
E(i) 




E(t) 



(3.19) 




which has a fixed point for (Q(t), cr(t), = (Q,0,0j, Q being a constant. Around this 
point the asymptotic behaviour is described by 



r Q(i) = 0, 

< a(t) = E(t), (3.20) 
±(t) = -V"(0)a(t) - 2QS(t) = -^ff(t) - 2QE(i). 
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Q decouples (and remains constant) and the only non-trivial remaining equation of motion 
is 



d 2 a(t) 4i? , da(t) 



(3.21] 



which can be integrated to give 



a(t) = C\ exp 



Q+ \ Q 2 - 4, J i 



+ C2 exp 



Q- \ Q 2 U J 1 



, (3.22a) 



Q + \ Q 



AR 



. dC 2 
+ — exp 



-Qt, 

(3.22b) 



where C\ and C2 are integration constants. 

For positive Q the solution converges to a = with or without oscillations depending 
on whether Q 2 5s AR/d, which clarifies the meaning of "large" and "small" Q(0) in Fig. 2. 

In terms of a(t) and 3>(t), this limit solution is 



a(t) 



t— >oo 







*(t) ~ $0 - Qt, 



(3.23) 



which is the conformal model in Eq. (2.1) plus a linear dilaton. 



Large a. In the large-cx limit, the solution goes to flat space. To study its behaviour let 
us consider the system 



In the large-a limit V(a) 



(Q(t) = -d{-V(a) + lQ(t) 2 ), 
\a(t) = -V'(a(t))-2a(t)Q(t). 

C and V'(a) — ► so it reduces to 

(Q(t) = -d{-C+lQ(t) 2 ), 
\a(t) = -2&(t)Q(t), 



(3.24) 



(3.25) 



which can be solved analytically: 




Q(t) = \l^ tanh (2VCdt + Ci) , 
a(t) = C3 + C2 arctan ^sinh C\ + cosh C\ tanh 
<D(t) = $ - J-log (cosh (tVCdt + cty + ^a(t) 




(3.26) 
(3.27) 
(3.28) 



where C\,Ci,Cj, are integration constants. As expected, for large t, both Q(t) and a(t) 
asymptotically converge to constant values and the background is flat with a linear dilaton. 
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3.3 The meaning of Q 

The parameter Q is interpreted as the background charge that drives the dilaton, in the 
regime where the latter is linear with time. This parameter dictates the behaviour of 
the solution around the fixed point in many ways. Apart from discriminating between 
exponentially decreasing and oscillating solutions, it also measures the consistency of the 
perturbative string approach. In fact positive values of Q correspond to a dilaton evolving 
towards large negative values, i.e. towards the perturbative regime. On the other hand, 
negative values of Q would lead to diverging solutions in which we would lose control over 
the approximations we made. It is worth to remark that in the non-linear dynamics Q(t) 
is not allowed to change its sign if we make the hypothesis of unicity of the solution. In 
fact changing sign would require to cross the Q = 0,a = const, line corresponding to the 
Euclidean conformal model. 

A final remark regards the consistency of the approximation for the dynamics one 
obtains from the RG-flow equation (2.9), corresponding to a Q — ► oo limit. The linearized 
system (3.20) provides a justification for such limit: in fact the time scale for Q(t) is 
comparably larger than the time scale for a(t) - to the point that the former decouples 
around the fixed point. For this reason it can be taken as a constant (fixed by the initial 
conditions) if we just concentrate on the evolution of the warp factor cr(t). 

3.4 The anisotropic deformations of the three-sphere 

We will now study the anisotropic deformation of the SU (2) model (see Sec. 2.2) when an 
extra time dimension is added to the system while the parameters 7^ (Eq. (2.13)) and the 
dilaton are allowed to depend on it. Explicitly, the full four-dimensional, time-dependent 
metric under consideration is 

3 

ds 2 = -dt 2 + J2^(t)J u J u = -dt 2 + fc((7 2 (t)cos 2 (/) + 7i(t) sin 2 0) d9 2 + 

+ (71 (t) cos 2 9 cos 2 (f) + 73 (t) sin 2 9 + 72 (t) cos 2 9 sin 2 <f>) d^ 2 + 
+ 73 (W 2 + 2 73 (t) sm6d4>d4> + 2 (71 (t) - 72 (t)) cos 9 cos 4> sin (f>d9di?j. (3.29) 

We introduce again the functions a v {t) as 

7 „(i) = e 2 ^W (3.30) 

and find the following equations of motion: 
2*(t) = £(Mi) + Mi) 2 ), 

v=l 

3 / 3 \ 

a^t) = R (l - e 4CT1 W + (e 2CT2 « - e 2 ^))^ f{ e^M - a^t) I -2*(i) + £ a u (t) I , 

(3.31) 
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plus permutations, where R = 6 /fc is the curvature of the undeformed three-sphere. Here, 
it is natural to define 



1 3 

Q{t) = -*(*) + 2 !>"(*) ( 3 - 32 ) 



and recast the system in the form 

3 

1 



■ &l{ t) = |(l - e 4CTl W + (e 2CT2 W - e 2 ^) 2 ] f[ e~ 2 ^ - 2& 1 (t)Q(t), 



(3.33) 



plus permutations. 

The analysis of Eqs. (3.33) can be performed in analogy with the system we obtained 
above, when considering only the overall breathing mode. We first observe that given the 
RG-flow equations 9 

^-/Wrt). (3-34) 
we obtain those for the time-dependent system as 

%^ = A,M*))-2W)^f (3.35) 

where Q(t) is the charge that naturally generalizes the one we found before. From these we 
conclude that the structure is again that of a dissipative system with the time-dependent 
friction coefficient given by the charge. We can further proceed by linearization around the 
577(2) wzw fixed point. This goes along the same lines as in the isotropic case and gives 

{* ) = \ s - (3.36) 

\a v {t) = -^a v {t)-2Qa v {t), u= 1,2,3. 

The charge Q(t) decouples and remains constant Q(t) = Q(0) = Q while the three modes 
become independent from each other and obey the same differential equation. This means 
that they oscillate in phase with each other and eventually converge to the fixed point 
corresponding to SU(2)k wzw plus a linear dilaton. 

All the considerations we made for the breathing mode then extend directly to this 
anisotropic generalization. The net effect of the time dependence is therefore to wash out 
the anisotropy, as was expected from the pure RG-flow analysis of Sec. 2.2. 



4. Cosmological interpretation 

The type of backgrounds we are studying are time-dependent and as such can be of cos- 
mological interest. For this reason, since there is a non-trivial dilaton, one should move 

9 Notice that f3 v (a) = M7)/2 7 „. 
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to the Einstein frame (as opposed to the string frame we have been using thus far). The 
metric takes the form 

~ a _ p -*(t)/2 s (A-,) 

i/MN — " 1/MN) K^- 1 ) 

and after a coordinate change 

r(t) = /V*(*')/ 4 di', (4.2) 
it reduces to the same warped product form as in Eq. (3.2): 
da 5 = 5 MN dx M dx N = -dr 2 + e Wt)-*{t)l* ^da^dx" = -dr 2 + a(r)g IJU ,dx li dx v (4.3) 

t=t(r) 



which has Hubble parameter: 

tf(T) = A log(a(r)) = e *M/4 U(t)-^ 



(4.4) 

t=t(r) 



Cosmologically interesting solutions are obtained when d = 3. In this case the Kalb- 
Ramond field i?[ 3 ] is proportional to the volume form on g. This implies that H 2 U oc g^ u 
and then Eq. (3.10) reduces to 

= 2A y> , (4.5) 

i.e. g^ u is the metric of an Einstein three-manifold. The simplest case is a three-sphere for 
which we get a typical example of FRW space-time of the kind studied in [36, 37, 38, 39]. 
As such it describes the time evolution of a homogeneous and isotropic 10 space-time, or 
more generally of a space-time with the target-space symmetries of the conformal theory 
in Eq. (2.1). 

Some intuition about the time evolution can be developed if we take the linearized 
system in Eq. (3.20). 

Linear Dilaton. In fact, as remarked above, if we consider the large t limit, the solution 
asymptotically approaches a linear dilaton background (which was already studied from 
this point of view in [14]): 

a(t) ► Q(t) = Q <f>(t) ~ -Qt, (4.6) 

t^oo 

hence one can easily verify that the metric in the Einstein frame is asymptotically 

cb? ^ _d T 2 + Q 2 T 2 g IIU dx^dx ,/ , (4.7) 
which corresponds to an expanding universe with Hubble parameter 

H(r) = 2 - (4.8) 



and curvature 

R + Q 2 d{d- 1 



3 We will discuss in this section the isotropic situation only. 
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Exponentially decreasing a. A similar result with a polynomial expansion is found if 
we consider an exponential decrease for cr(t), or better for c(t) (in the linear limit c(t) — 1 
obeys the same equations as a(t)) i.e. if we are in the large-Q regime. After a redefinition 
of the variables we can write: 

c(t) = e~ at + 1. (4.10) 

It is easy to check that in general 11 

r(t) = Jc(t)- d ^ 4 f^ dt 'dt, (4.11) 
and in this linearized approximation this becomes 

r(t) = J (e-' + lj-^e^dt. (4.12) 
Although this integral can be solved analytically: 

(4.13) 

where u = e* and 2 -Pi is an hypergeometric function 12 . It is better to consider the asymp- 
totic behaviour. For u — > 00 one finds that r(u) and the warp factor a(ii) go as: 

T (u) ~ i^ /4 , o(«) ~ ^ /2 , (4.16) 

t^oo (J r^oo 

and consistently with the results above for the linear dilaton case (which is precisely the 
large- u limit): 



r^oo 

similarly for small u: 



a(r) ~ t 5 (4.17) 

T^OO 



r(«) ~ ^- (d+40)/16 , a(«) ~ „-/^/» (4.18) 

and then 

a(r) ~ r 4+8(4-Q)/(d+4Q) . (4 . 19) 

Since we have a power-law expansion at both limits, the Hubble parameter goes as H{t) oc 
1/t and more precisely is asymptotically given by: 

H(r) ~ 2 tf(r) ~ - fl + ^iV (4.20) 



11 Since c(t) > by construction, the relation r = r(t) is always invertible. 
12 The hypergeometric function 2 f\ is defined as follows: 



where (a) k is the Pochhammer symbol 



2Fl (aAc,z) = ±^^i, (4.14) 



r(q + fc) 

(a) * = i^r- (4 - 15) 



- 19 - 



Note that this behaviour precisely measures the effect of a finite value for Q and in 
fact for Q — > oo we recover again H(t) ~ 2/r. To summarize, we get again a polynomially 
expanding universe (a so-called big-bang solution). 

Damped oscillation regime. The small-Q regime is more difficult to be studied ana- 
lytically. Nevertheless we can make some qualitative comments on the overall behaviour 
of the solution. 

First of all we can redefine the variables so to recast Eqs. (3.22) in the form 



a{t) = Ae~ Qt cos(Lot) 



Ad 



<f>(t) = $ + ^e~ Qt cos(wt) - Qt 



(4.21a) 
(4.21b) 



where A is a constant and uj 2 = \Q 2 — 4R/d\. The warp factor and the Hubble parameter 



then read 



a(r) 



H(t) 



a 2«r(t)-*(t)/2 



a *(t)/4 



t=t(r) 



exp 



Ae-& cos(wt) ( 2 





= exp 




t=t(r) 



4 



Ad 



-Qt 



$0 Qt 
2 2 

cos(wt) - 



-Ae _( ^ (Qcos(u;t) +cjsin(^)) ^2 - ^ + ^ 



t=t(r) 

Qt' 

4 



t=t(T) 



(4.22) 



(4.23) 



where r(t) is defined by 



r(t)= r e -*(*')/ 4 di'= f 



exp 



- T -— e 0* cosM) + — 



di'. 



(4.24) 



The Einstein frame time variable r(t) is given by the integral of a positive function of t. 
This means that it is monotonically increasing and the relation r = r(t) is invertible. In 
turn the qualitative behaviour of the functions a = a(r) and H = H{t) can be understood 
by studying their behaviours as functions of the string frame time t. In particular H(t) 
is an oscillating function and has a zero for each zero of the function (^2a(t) — &(t)/2^j. 
The maximum number TV of such zeros can be estimated in the following way. If we write 
explicitly the equation 2a — 3>/2 = we obtain 

d 



Q = 2Ae~ Qt (Q cos(u;t) + uj sm(ojt)) 2 



(4.25) 



so, we are looking for the intersections of the line w(t) = Q with the function w{t) 
2Ae-& (Qcos(l^) + usm(ut)) (2 - f). The latter is bounded by 



w(t) = 2Ae~ Qt ^Q 2 + 



and then intersections are only possible for t < t where t is given by 

_ 1 A^Q 2 + u 2 |8 - d\ 



Q 



log 



2Q 



(4.26) 



(4.27) 
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Figure 3: Typical behaviour for the warping factor and Hubble parameter in the small-Q regime. 
Numerical solution with parameters A = l,Q = 1.2, ui = 10, d = 3. 



In the (0, t) interval, the (Q cos(cot) + to sin(u;i)) term can oscillate cot/ (2ir) times and then 
the number of solutions for Eq. (4.25) is at most 



n < N 



LOt 




TT 





to , AJQ 2 + u 2 18 - d\ 

log — = — 



(4.28) 



The overall behaviour for the warp factor a(r) is of that of a power-law expansion super- 
imposed with n < N oscillations. The Hubble parameter after n oscillations goes to zero 
for large r as H(t) ~ 2/r. Typical behaviours for a(r) and H{t) are shown in Fig. 3. 



5. Summary 

Finding viable time-dependent solutions with cosmological interest is one of the open prob- 
lems in string theory. In this note we have presented a family of such solutions, derived as 
deformations of a conformal model, leading to a generalization of frw universe undergoing 
a big bang-like expansion superposed, in some regimes, with damped oscillations. In these 
backgrounds the time dependence comes from a Liouville-type mechanism in which flow- 
ing from the conformal point is seen as the effect of a time-dependent field. The RG-flow 
equations on constant-time slices then become an approximation for the Weyl-invariance 
equations and the off-shell dynamics is interpreted as the on-shell dynamics with one extra 
dimension. As already stressed, this behaviour is not generic and the identification of the 
RG-flow with time evolution is not valid in any regime. 

The constraints imposed by the rational nature of the underlying compact model es- 
sentially fix the type of deformations one can consider and naturally lead to an isotropic 
FRW background whose behaviour is due to the presence of an axion (dual to the Kalb- 
Ramond field) and an oscillating dilaton. Anisotropic extensions of FRW solution are also 
considered but time evolution smoothes the initial anisotropy and the system relaxes to a 
symmetric phase. 

A common feature of all situations considered here is that the RG equations at lowest 
order in a' always describe a Ricci flow with torsion. In the case of the three-sphere with 
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anisotropic perturbation, the Ricci-flow equations are identical to the Halphen system, and 
are thereby integrable by using modular forms. More general perturbations exist, which 
are compatible with the mini-superspace approximation (see [20] ) but a systematic analysis 
of all those is beyond the scope of this work. 

It stems from the overall analysis performed in the present framework that in a quite 
universal manner, sooner or later, time evolution is identified with an RG flow, which is 
a Ricci flow. Time becomes thus a two-dimensional scale, the dilaton being the interplay 
between target space and worldsheet. Thurston's geometrization conjecture is then at 
work. It tells us that the string target space will universally converge towards a collection 
of Einstein spaces, whenever these are available, enhancing therefore the symmetry. This 
behaviour is independent of the initial conditions. 

Clearly, this reasoning is valid under the assumption that the initial space is locally 
homogeneous. Exact time-dependent string backgrounds with non-homogeneous spaces 
exist [40, 41], and there time evolution does not lead to Einstein spaces. Consistently with 
the above picture, in these cases, time evolution never identifies with an RG flow because 
it is rather related to a marginal deformation. 

The latter observations are puzzling additions to the already mysterious origin of time 
in string theory, along the lines of thought of A. Polyakov. 
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A. The Ricci tensor for the deformed group manifold 

In order to calculate R a B, we need to use some notions in group manifold geometry. Let 
{ a } be a set of one-forms on a manifold Ai satisfying the commutation relations 

[& ) M = f a fh8 a (A.1) 

as it is the case when a are the Maurer-Cartan one-forms for a Lie algebra and f a ^ 
the corresponding structure constants. We wish to study the geometry of the Riemann 
manifold M endowed with the metric 

9 = 9aie a ®^. (A.2) 
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In general such a metric has a symmetry G xG' where G is the group corresponding to the 
structure constants /°g 7 and G' C G. The maximally symmetric case, in which G' = G is 
obtained when g is G-invariant, i.e. when it satisfies 

f%9aS + f\g a p = 0. (A.3) 

For compact groups this condition is fulfilled by the Killing metric. 

The connection one-forms io a ^ are uniquely determined by the compatibility condition 
and the vanishing of the torsion. Respectively: 



dg a/ 3 - uj 1 ^^ - uP^a = 0, 
d§ a + u a p A ¥ = T a = 0. 

As it is shown in [42], if g a p is constant, the solution to the system can be put in the 

ufp = -D°^ 

where D a ^ = 1/2 f a ^ — K a ^ and K a p is a tensor (symmetric in the lower indices) 
by: 

The Riemann tensor and the corresponding Ricci tensor now read: 
n 13-yS — u /3kJ 7<5 + u k^ U 05 ~ U k& u /3 7 

R/3S = D a /3K f K aS - D a KS D K p a . 
B. The Jacobi functions 

We take the standard conventions for the Jacobi functions: 

pel 
pel 
pel 

where q = exp2iirz. With these conventions 

A-** 1 A — ^ 

and 

2rf = 2 3 4 

with 

00 

V (z) = qh. H (1 - q n ) 
n=l 



A.4 
A.5 

form 

A.6 

iven 

A.7 

A.8 
A.9 



B.l 



B.2 



B.3 



B.4 



B.5 



B.6 
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the Dedekind function. The non-holomorphic function 

12^ 3 
in az 7rlm z 

is modular covariant of degree 2. The holomorphic part of the latter, which is not modular- 
covariant, is the function y (divided by iir) in (2.31). 
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